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A first-principle theoretical approach to study the process of radiative electron attachment is de¬ 
veloped and applied to the negative molecular ions CN“, C 4 H“, and C 2 H“. Among these anions, 
the first two have already been observed in the interstellar space. Cross sections and rate coef¬ 
ficients for formation of these ions by radiative electron attachment to the corresponding neutral 
radicals are calculated. For completeness of the theoretical approach, two pathways for the process 
have been considered: (i) A direct pathway, in which the electron in collision with the molecule 
spontaneously emits a photon and forms a negative ion in one of the lowest vibrational levels, and 
(ii) an indirect, or two-step pathway, in which the electron is initially captured through non-Born- 
Oppenheimer coupling into a vibrationally resonant excited state of the anion, which then stabilizes 
by radiative decay. We develop a general model to describe the second pathway and show that 
its contribution to the formation of cosmic anions is small in comparison to the direct mechanism. 
The obtained rate coefficients at 30 K are 7 x 10~^®cm®/s for CN~, 7 x 10“^^cm®/s for C 2 H“, and 
2 X 10“^®cm®/s for C 4 H“. These rates weakly depend on temperature between lOK and 100 K. The 
validity of our calculations is verified by comparing the present theoretical results with data from 
recent photodetachment experiments. 

PACS numbers: 52.20.Fs, 32.80.-t, 32.80.Fb, 33.80.-b, 34.80.Lx 


I. INTRODUCTION 

The present theoretical study of radiative electron at¬ 
tachment (REA) to neutral molecules is motivated by the 
recent discoveries of molecular anions in the interstellar 
medium (ISM). Six anions have been detected so far in 
the ISM: CeH- C 4 H- 0, CgR- 0, CgN" 0, 

C 5 N- 00, andCN- 0 . The possibility for atomic an¬ 
ions, such as H~, to be formed in the ISM by REA was 
first suggested by McDowell in 1961. Later, Dal- 
garno and McCray [ll| have discussed the role of nega¬ 
tive atomic ions in the formation of neutral molecules in 
the ISM. The formation of molecular anions in the ISM 
by REA has been proposed by Herb st who has also 
developed a theoretical approach 0, Il 2 || to evaluate rate 
coefficients for REA. More than twenty years after his 
prediction, negative molecular ions were indeed detected 
in the ISM. 

The theoretical approach proposed by Herbst, the 
phase-space theory (PST), has been used in a number of 
studies 01 IM 1 to calculate the REA rate coefficients 
and to model formation of anions in the ISM. The ap¬ 
proach has been successful in interpreting the observed 
column density of CsH”, CeH”, CsN” and CsN” ions, 
while the agreement with observations is not as good for 
the C 4 H“ and CN“ ions. 

PST relies on several assumptions. First, it considers 
REA as a two-step process, schematically represented in 
Fig. [T]for the case of CN/CN“. As a first step, the inci¬ 
dent electron is captured by a neutral target molecule 
M into an electronic state of M~ through non-Born- 


Oppenheimer coupling, thus forming a vibrational res¬ 
onance of M~ in a highly-excited vibrational level. The 
resonance can decay back to the M+e~ electronic contin¬ 
uum spectrum through the same non-Born-Oppenheimer 
coupling, i.e. the electron can autodetach. Another pos¬ 
sibility for the resonance, considered as the second PST 
step, is to emit a photon, thus stabilizing the M~ system. 
In PST, it is assumed that the probability of the first step 
of the process is unity and, therefore, the cross section 
for the electron capture is approximated by the unitary 
limit formula for s-wave scattering 0 cTc = where 

k is the wave number of the incident electron with energy 
Eel = {Kk)^/{2me), and rrie is the electron mass. In the 
second PST step, stabilization of the M~ resonance is 
represented as an emission of a photon by a set of har¬ 
monic oscillators of the molecule in the normal mode ap¬ 
proximation of molecular vibrations. A larger number of 
available vibrational modes decreases the probability of 
autodetachment, and thus allows for a larger probability 
in the second PST step. Therefore, in PST, the over¬ 
all two-step REA cross section apsT grows rapidly with 
the number of atoms in the molecule and approaches the 
unitary limit upsr For example, the REA rate co¬ 

efficient for formation of CN“ calculated by PST is about 
10 “^^cm^/s [a and much larger, 3 x 10 ^cm^/s, for 
CeH- 0 at lOK. For CN the theoretical value is smaller 
by several orders of magnitude than the one needed to 
explain the [CN“]/[CN] abundance ratio obtained from 
the astrophysical observations. 

The PST assumption about the unitary probability 
of the first step of the process can hardly be justified 
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for small molecular ions. But it is possible to apply a 
quantum-mechanical approach based on first principles 
for the both steps of the REA mechanism, suggested by 
Herbst. Considering the process quantum-mechanically, 
another mechanism for REA is also possible: The elec¬ 
tron incident on a neutral molecule emits a photon and 
becomes bound without an intermediate step of capture 
into a vibrationally excited state of the ground electronic 
state of the negative ion. To distinguish the two mech¬ 
anisms, we call the first one as indirect (IREA) and the 
second mechanism as direct REA (DREA). The total 
cross section of the REA is the sum of DREA and IREA 
cross sections. 

Recently, we have developed a fully-quantum theoret¬ 
ical approach for DREA based on first principles only 
[l^ . The approach considers the radiative electron at¬ 
tachment of a continuum electron through spontaneous 
decay to the anion ground state (see Fig. [T]) and does not 
include an intermediate vibronic state of M~ populated 
through the non-Born-Oppenheimer coupling, as repre¬ 
sented by the two-step REA mechanism in PST. Hence, 
the approach relies on wave functions of the continuum 
spectrum of the M + e~ system and transition dipole mo¬ 
ments to the bound electronic state of M“, calculated ab 
initio. Using the developed approach, we have calculated 
the REA cross section and rate coefficient for the for¬ 
mation of the cyanide ion, CN“. Our results confirmed 
the previous assessment that the REA rate coefficient for 
CN“ is too small, 8 x 10“^®cm^/s at 10 K [l^, to explain 
the CN“ abundance observed in the ISM. 

In the present study, we extend the theoretical ap¬ 
proach to larger molecules. In order to assess the ap¬ 
proximations employed in the PST approach, we also de¬ 
velop a quantum mechanical approach for IREA. In this 
study, we calculate explicitly the probability of electron 
capture during the first step of IREA process using ab 
initio methods and determine the overall cross section of 
the IREA process, which might be compared with the 
PST results. 

Although there is no experimental data on the REA 
process for carbon chain molecules, using a similar the¬ 
oretical approach, we determine in this study cross sec¬ 
tions for the inverse process to REA, namely photode¬ 
tachment (PD), and compare with available data from 
recent photodetachment experiments [H, This al¬ 
lows us to verify the validity of our results. 

In the next section, we present our theoretical ap¬ 
proach to study DREA and apply it to the CN“, C 2 H“, 
and C 4 H“ ions. Section imi is devoted to the compari¬ 
son of the results obtained in this study with data from 
photodetachment experiments. Section IIVI presents the 
theoretical approach of IREA for the case of CN“. In 
the Section El we develop a model of IREA for larger 
molecules and we compare its results with those of PST. 
Finally, Section ED summarizes the important findings of 
the study. 



FIG. 1: Schematic representation of two mechanisms for REA 
of an electron to the CN radical: (i) DREA: The electron 
incident on CN in its ground vibronic state spontaneously 
emits a photon of energy (green arrow) and forms the 
CN~ ground state. The photon energy is equal to the dif¬ 
ference between the initial total energy of the system Etot 
(horizontal dashed line at near energy of the v~ = 19 state) 
and the energy of CN“ ground state, (ii) IREA: As a first 
step of the process, the incident electron is captured via non- 
Born-Oppenheimer coupling into the ground electronic state 
of CN“ without emitting a photon. Because the total energy 
Etot is conserved, the vibrational level v~ of CN“ in this first 
step of the process is highly excited: For low incident energies 
of the electron, it corresponds to v~ = 18 or 19. After the 
first step, the electron captured in an excited vibrational level 
can either autodetach or stabilize by photon emission. The 
photon changes the rotational and vibrational states of the 
CN“ molecule. The largest probability for the second step 
corresponds to a change of v~ by one quantum. The dashed 
red line shows the derivative of the vibrational ground state 
wave function of CN, which enters in the calculation of the 
non-Born-Oppenheimer coupling in Eq. mi)- 


II. DIRECT MECHANISM OF RADIATIVE 
ELECTRON ATTACHMENT TO CN, CaH, AND 

C4H 

The cross section for DREA of an electron to a neu¬ 
tral linear molecule M, such as Cn,H (n = 2,4) or C^N 
(to = 1,3), initially in its electronic ground state with 
the vibrational level v and energy Et, was given in [l8l | 
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FIG. 2: Matrix elements of the and components of 
the transition dipole moment between the C 2 H“ electronic 
ground state and the e“+C 2 H system for several partial waves 
as a function of the incident electron energy, calculated at the 
equilibrium geometry of C 2 H~. 


and is expressed as 

Itt 


( 1 ) 


where Vf is the vibrational state of the ion M~ with to¬ 
tal energy Ef formed after DREA; oj is the frequency of 
the emitted photon, huj = Ei + E^i — Ef. The quan¬ 
tities are matrix elements of the components 

TT = —1,0,-|-1 of the dipole moment operator between 
the initial iPp/A (M-bincident electron) and the final iP/ 
electronic state, integrated over the initial Xvio) and final 
Xvfio) vibrational wave functions of M and M“, respec¬ 
tively, and q denotes collectively all internuclear degrees 
of freedom. Their values are given by 

Xvj{^{'i’f\d^\'i’ri\)Xv{^dq, ( 2 ) 

(see details in Ref. My, with I the electronic partial 
wave angular momentum and A its projection on a spe¬ 
cific axis in the molecular frame. The matrix elements 
(d'/ld^l'I'riA) of the dipole moment operator are given 
by the integral 


N 

(T/ld^ld'ru) = “X! / ,'^N)erk-„ 

k=l 

xd'riA(ri, • • • ,rAr)d^ri • • • d^CAT, (3) 


where the function lb f represents the fV-electron final 
state of the negative ion (ri,...,rAr are the coordinates 


of the electrons), dtpiA is the electronic continuum state 
representing the scattering electron with I and A angular 
quantum numbers, and T labels the initial neutral elec¬ 
tronic target state with — 1 electrons. Finally, is 
one of the three cyclic components (tt = 0, ±1) of the 
coordinate of the electron 


_ J Zfc, if tt = 0 

\ T{xk±iyk)/V2, a TT = ±1. 

The calculations of the electronic wave functions and 
transition dipole moments (TDMs) of Eq. ([3]) are per¬ 
formed using the complex Kohn variational method, ex¬ 
tensively described in past studies [U . 

The DREA cross section for the formation of CN“, 
starting from the ground vibrational level u = 0 of CN, 
was calculated using Eq. O 0. The vibrational in¬ 
tegral of Eq. @ was computed explicitly from the 
geometry-dependent matrix elements ('f'/ld^l'I'riA) ob¬ 
tained in the complex Kohn calculations. These matrix 
elements, as a function of the electron energy and of the 
internuclear distance, were respectively shown in Fig. 3 
and Fig. 4 of Ref. [II. 

The calculation of the TDMs is computationall y in ten- 
sive, especially for polyatomic molecules. In Ref. [I^, we 
found that the TDMs weakly depend on the geometry of 
the molecule near its equilibrium position. Therefore, 
it seems reasonable to use the Franck-Condon approx¬ 
imation and simplify the calculation of the vibrational 
integral in Eq. (ED by using the value of (4>/|d^|4>riA) at 
a fixed molecular geometry, e. g. the equilibrium posi¬ 
tion of the negative ion or the equilibrium of the neutral 
molecule, which are close to each other. The TDMs then 



energy (eV) 

FIG. 3: Same as Fig. [5] for the e“-|-G 4 H system. 
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take the simple form 

~ (^/M^l^nA)eo J xlfi^Xvi^dq, (5) 

where the subscript Qq refers to the equilibrium geome¬ 
try of the negative ion. For molecules such as CnH and 
CmN, for which the potential energy surfaces of the ini¬ 
tial electronic state of the target and the final state of the 
negative ion are quite similar in shape near the equilib¬ 
rium positions of the ion and the neutral molecule [^ . 
the Franck-Condon integral in Eq. ([S]) is the largest for 
transitions with Vf = v, for which its value is close to 
unity. For instance, the Franck-Condon integral in Eq. 
® is about 0.90 for C 2 H/C 2 H“ and 0.87 for C 4 H/C 4 H“ 
|23l |. For transitions to other vibrational levels, the inte¬ 
gral is significantly smaller. Therefore, the DREA cross 
section is well approximated by 

^ ^ ^ : ( 6 ) 

Itt 

where the transition dipole moment is evaluated at the 
energy of the M + e~ system. The transition dipole mo¬ 
ments strongly depend on energy, especially for non-zero 
partial waves, 1 > 0. At low energies, their behavior is de¬ 
scribed quite well by the Wigner threshold law [2^ . The 
energy dependence of |('k/|d^|'I'riA)eol for C 2 H/C 2 H“ 
and C 4 H/C 4 H“ is shown in Eigs. [2] and [3 

For C 2 H and C 4 H, the DREA cross sections were calcu¬ 
lated using the approximate formula of Eq. ®, with the 
transition dipole moments determined only at the equi¬ 
librium geometry Qq of the negative ion for several ener¬ 
gies of the incident electron. Eor comparison, the DREA 
cross section for CN“ has been calculated using the di¬ 
rect integration over internuclear distances, Eqs. o and 
as well as using Eq. (jS]). Figure 0] shows the DREA 
cross sections obtained for CN“, C 2 H“, and C 4 H“. As 
can be seen in the figure, the CN“ cross sections ob¬ 
tained in the two ways are almost identical. Therefore, 
Eq. ([H]) gives a very good approximation of Eqs. (ED 
and ([2]) for CN . For C 2 H/C 2 H and C 4 H/C 4 H the 
approximation is likely less accurate because the vibra¬ 
tional functions of the neutral molecule and the ion are 
not as similar to each other as for CN“. 

The obtained DREA cross sections have been used to 
determine the thermal rate coefficients, which are shown 
in Eig. [5l The rates for the formations of these three an¬ 
ions via DREA are too small to explain the astrophysical 
observations. 


III. COMPARISON WITH THE RESULTS OF 
PHOTODETACHMENT EXPERIMENTS 

There is no experimental data on radiative electron 
attachment to the CN, C 2 H, and C 4 H molecules. How¬ 
ever, the calculated transition dipole moments can be 
used to determine the photodetachment cross sections. 



FIG. 4: The cross sections for formation of CN“, C 2 H“, and 
C 4 H“ by DREA with the ground vibrational level a = 0 of the 
neutral molecule as the initial state. The CN“ cross section is 
calculated in two different ways: The dotted line corresponds 
to the calculation using Eq. the cross section shown with 
the solid line, ”vib FT” (vibrational frame transformation), 
is obtained evaluating the vibrational integral of Eq. m ex¬ 
plicitly. 



FIG. 5: The REA thermally-averaged rate coefficients. 


for which experimental data on absolute values of pho¬ 
todetachment cross sections have recently been obtained 
for CN“ [ 1 ^, C 2 H“ [13 and C 4 H“ [ 2 ^ anions. Here, we 
discuss the interpretation of the photodetachment experi¬ 
ments only briefly. A detailed and more elaborated study 
can be found in Ref. [^ . 

The photodetachment cross section is given by (see, 




















5 



FIG. 6 : Theoretical (lines) and experimental (symbols with 
uncertainty bars) photodetachment cross sections for CN~, 
C 2 H~, and C 4 H~. The triangle is the experimental result for 
CN- [ii], diamonds - C 2 H- [13] , squares - C 4 H- [ 23 |. 


for example, Eq. (2.202) of Ref. phi ]) 


crpD 


3c 


\di,Tr,E\^ , 
Ztt 


( 7 ) 


where u) is the photodetached photon frequency and 
di^TT,E is the dipole moment between the initial bound 
state and the final continuum state with energy E and 
partial wave I of the photodetached electron. The radial 
part of the wave function (j)E{f) of the initial electronic- 
continuum state used to calculate the transition dipole 
moment (Ii^t^^e in Eq. 0 is energy-normalized. The con¬ 
tinuum functions used in our calculation (see Eq. (|3|)) 
are normalized as 


, . /2me 

^E[r) = 


( 8 ) 


Thus, the PD cross section can be written as 


SrUeTru) 


apD = 




Itz 


-^7r,r/—TT 


(9) 


Eor convenience of comparison with the results of the 
DREA calculations, the photon frequency can be ex¬ 
pressed in terms of the electron affinity Eea and the en¬ 
ergy Ef .1 of the incident (in DREA) or emitted (in PD) 
electron huj = E^i -\- Eea- This formula assumes that 
both initial and final vibrational levels are not excited, 
V = v~ = 0, and that the zero-point energies for the 
neutral molecule and the ion are the same. The exper¬ 
imentally measured affinities are 3.862 eV for C N f^ . 
2.969 eV for C 2 H [Hil, and 3.558 eV for C 4 H 
We discuss now the accuracy of the present calcula¬ 
tions. The main source of uncertainty for the REA and 
PD cross sections for all three molecules considered here 
is the quality of electronic continuum and bound state 


wave functions. The quality of the electronic bound state 
wave function could be assessed in part by comparing the 
obtained theoretical affinity with the experimental one. 
The affinities corresponding to wave functions used in 
the present calculations, are 3.8 eV, 2.2 eV, and 3.0 eV 
for CN“, C 2 H“, and C 4 H“ respectively. Therefore, the 
agreement is about 1% for CN“ and much poorer, about 
30%, for C 2 H“ and C 4 H“. The quality of electronic con¬ 
tinuum wave functions is more difficult to assess. Based 
on our previous experience with the electron-scattering 
calculations, we assume here that an additional uncer¬ 
tainty in the calculated transition dipole moments due 
to the quality of continuum wave functions is about 10 % 
for CN“, C 2 H“, and C 4 H“. These considerations give an 
estimated uncertainty in the REA and PD cross sections 
of about 20% for CN“ and 40% for C 2 H“ and C 4 H“. 
Eor C 2 H“ and C 4 H“, there are additional sources of 
uncertainty; the neglected geometry dependence of the 
transition dipole moments and the neglected role of rovi- 
brational Feshbach resonances, which could be present at 
low collisional energies. 

Figure | 6 ] shows the PD cross sections calculated for 
CN“, C 2 H“, and C 4 H“ using Eq. (0 and compares 
them with the available experimental data, which were 
estimated in Ref. to have about 25 % of uncertainty. 
The agreement is very good for C 2 H“, especially for ex¬ 
perimental data points near the photodetachment thresh¬ 
old. The agreement for C 4 H“ is also good for the lowest 
energy point, but is about a factor 2-3 lower than the ex¬ 
perimental value for the second energy point. Similarly, 
the only experimental data point for CN“ measured at 
4.65 eV gives a PD cross section twice larger than the 
theoretical value. The reason for this discrepancy is not 
clear. Overall, the agreement of the theoretical and ex¬ 
perimental results is sufficient to conclude that the theory 
is reliable for calculations of photodetachment cross sec¬ 
tions and, therefore, also for DREA cross sections. Thus, 
the results of photodetachment experiments validate the 
present theoretical approach for the DREA process and 
the results obtained using the approach. In particular, 
the present results and the PD experiments [i3,l23 
C 4 H“ and CN“ suggest that the observed abundance of 
these two ions in the ISM can hardly be explained by the 
DREA mechanism. 


IV. THE CONTRIBUTION OF THE INDIRECT 
PROCESS TO THE TOTAL REA CROSS 
SECTION: THE CN EXAMPLE 


We now consider the process of radiative attachment 
mediated by the non-adiabatic couplings, through the 
IREA mechanism discussed in the introduction. In this 
section, we consider the IREA process for CN, for which 
the non-adiabatic couplings are evaluated numerically. 
We extend the approach to larger molecules in the next 
section. 

The probability per unit time of a transition from the 
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initial vibronic state |i) of the e“+CN system described 
above into the final state of CN“ being in a vibrationally 
excited level n/ ~ 19 is given by the Fermi’s golden rule 


P = 


2tt 


(/|AK) p{E,), 


( 10 ) 


where p{Ec) is the density of final states after the electron 
capture p{Ec) ~ 1/AErv, with AErv the energy splitting 
between rovibrational states, and A is the operator of 
non-adiabatic coupling. Denoting R the CN internuclear 
distance, the matrix element A fi between initial and final 
vibronic states is 

= ^ I x:,{R)Hrix{R)^Xv{R)dR, ( 11 ) 



where pcN is the reduced mass of CN and Af^rixiR) is 
the following electronic matrix element 

Af,rix{R) = j «'}(f,i?)^«'r/A(r,i?)dr. (12) 

In the above equation, r denotes collectively all electronic 
coordinates. Note that we neglected the second term of 
the non-adiabatic couplings, which is expressed with the 
second derivative of the electronic wave function with 
respect to the internuclear distance. 

Calculating A/i, we did not account for the integral 
over rotational coordinates. The latter integral is of the 
order of unity or smaller. A change of rotational quantum 
number j during the process of the non-adiabatic electron 
attachment is unlikely for CN“ at low collision energies. 
This is because the values of the non-adiabatic couplings 
are the largest for s-wave scattering and much smaller 
for higher partial waves at low energies. Thus, we use 
the vibrational splitting in the formula for the density of 
states p{Ec) ^ 1/AEy. 

Numerical evaluation of the electronic matrix ele¬ 
ment Af^Tix{R) of the non-adiabatic couplings was per¬ 
formed at several internuclear distances. Typically, non- 
adiabatic couplings can be calculated using standard ah 
initio programs, but such calculations are limited to cou¬ 
plings between electronic bound states, whereas in the 
present case, the initial electronic state belongs to the 
electronic continuum spectrum of the molecule. One pos¬ 
sible way to calculate the couplings with a continuum 
state using standard ab initio programs is to include very 
diffuse orbitals in the bound-state calculations in order 
to cover the region of large electronic radial distance. 
As more diffuse functions are added to the basis set, a 
better description of the asymptotic region is achieved, 
which leads to the appearance of "box-state” like wave 
functions with positive asymptotic energy. Such states 
resemble a scattering state of the electron if enough dif¬ 
fuse functions are added to the basis set. The box states 
should then be normalized appropriately to represent a 
continuum state. The latter approach should provide a 
reasonable approximation of the non-adiabatic couplings 
because no strong boundary conditions are imposed on 


FIG. 7: The matrix element of the non-adiabatic coupling 
operator of Eq. m evaluated for two different electronic 
"quasi-continuum” states with energies -^10 and ~60 meV 
above the electronic threshold. The continuum state of the 
incident electron is normalized in the same way as in Ref. [l^ , 
where the current density of electrons in the incident plane 
wave is equal to the electron velocity (see Eqs. (10), (A2), 
and (A19) in [H). 


the wave function. Performing numerical calculations, 
we have verified convergence of results (the non-adiabatic 
couplings) with respect to the number of added diffuse 
functions. 

Since the calculations are performed in a finite volume 
represented by the space spanned by our diffuse basis set, 
the calculated electronic states of CN-|-e“ just above the 
CN threshold energy appear as a "quasi-continuum” of 
states T'Jj with discretized energies. These states are la¬ 
beled (in order of increasing energy above the electronic 
threshold) with index n and with the dominant angular 
momentum 1. The states should be rescaled to be energy 
normalized by multiplying by a factor AE^^ ' , where 
AEei is the energy difference between ” quasi-continuum” 
states. Because AEei is changing from one state to an¬ 
other, it is taken here as the average between two en¬ 
ergy differences for three neighboring ” quasi-continuum” 
states. The calculated values of Af^rix{R) are plotted 
in Fig. [7] for two "quasi-continuum” s-wave scattering 
states. These states, corresponding to n = 2 and 3, and 
have respectively ^10 and ^60 meV asymptotic electron 
energy above the electronic threshold. Notice that the 
non-adiabatic couplings for the two states are similar 
and decrease with the internuclear distance. Further¬ 
more, states with higher energies have larger couplings, 
which is expected from the Wigner threshold law because 
such couplings should increase approximately as \/Eei. 
Because the low-energy non-adiabatic couplings associ¬ 
ated with p-wave scattering states increase only as 
their values at low energy were found to be negligible 
in comparison with the couplings from s-waves scatter¬ 
ing states. Finally, although the non-adiabatic couplings 
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were obtained from an approximative treatment, their 
values should represent an estimation accurate enough 
for the purpose of this study. 

The electron capture cross section is then obtained by 
dividing the probability by the current density jcd in the 
incident wave. If the incident wave is normalized as in 
Eqs. (A2) and (A19) of Ref. the current density is 
then simply equal to the velocity of the incident electron 
jcd = Vcd- It gives the following estimation for the cross 
section 


CTc 


P _ TT I 

jcd Eg/ APu 


X*vf (^)A/.rzA {R)dR 


(13) 


In the above equation, we assumed that the ratio of the 
number of final rotational levels of the formed anion CN“ 
to the number of initial rotational levels is approximately 
unity. Moreover, we consider that the value of Af^ri\{R) 
is almost constant in the energy range under study. This 
represents a relatively good approximation, at least be¬ 
tween 10-100 meV, on the view of the weak energy de¬ 
pendence of Af^rix{R), as seen in Fig. [71 

We have calculated the vibrational integral numerically 
with the final vibrational level of CN“ (u/ = 19), which 
has an energy close to the energy of the initial vibrational 
level of CN {vi =0) with a negligible asymptotic energy 
of the incident electron. We obtained the value 


/ O 

Xi 9 idi)Af^rix{R)-^Xo{R)dR = -1.5 x 10"'’a.u., 

(14) 

expressed in a.u. (atomic units). The vibrational energy 
spacing AEy is about 0.01 hartree (H) and the reduced 
mass ficN is 1.1 x 10"* rUe (a.u.). Therefore, the cross 
section to capture an electron into a vibrational level is 
approximately 




6 X 10-16 , 


(15) 


The estimate (Ha is an upper bound for the IREA 
cross section because of the autodetachment process: The 
actual cross section is reduced because the formed CN” 
ion in the excited vibrational level can decay back to 
the CN molecule and a free electron. The overall cross 
section of the indirect REA is then 

(JIREA = -^CTc < CTc , (16) 

tot 


where Ttot is the total width for the decay of the unsta¬ 
ble vibrational state of CN” formed during the collision. 
For the case of CN/CN”, the total width is the sum 
Ttot = Esp -I- Tad of the widths towards autodetachment 
Tad and towards spontaneous emission into all possible 
vibrational levels. 

The rate coefficient Tad can be calculated using the 
Fermi’s golden rule and the same matrix element of the 


non-adiabatic coupling of Eq. (HU). The density of final 
states is calculated differently since it will correspond to 
an outgoing electronic state. In this case, the density 
of states is unity because the radial part of the wave 
function is energy-normalized. Therefore, the probability 
of autodetachment per unit time is roughly 

Stt 

Tad ^ 

if, at the end of the process, the CN molecule is again in 
the ground vibrational level. 

The rate T^p of spontaneous emission can also be es¬ 
timated using the standard formula (see Eq. (2.189) of 
Ref. HI) 



TT 


where dlj’ ^ is the vibrational matrix element of the 
permanent dipole moment of the anion (R) for vibra¬ 
tional transition v~ ^ v' : 

j xl'-{R)d.{R)xv-{R)dR. (19) 


In the above equation v'~ is the vibrational level of CN” 
after emission of a photon, which should be smaller than 
v~ in order to stabilize the anion. The vibrational matrix 
element are the largest for Av~ = 1. For such a tran¬ 
sition, the vibrational dipole matrix element is roughly 
equal to the derivative of d-j^iR) with respect to R. The 
order of magnitude of the derivative is about unity in a.u. 
and, therefore, the vibrational matrix element of the per¬ 
manent dipole moment is about coq. For transitions with 
Av~ > 1, the matrix elements are significantly smaller. 
Thus, we only account for the v' = v~ — 1 transition, 
for which Wgp = AEy/h, and obtain T^jp ^ 1.3 x 10“^^ 
H/fi = 0.54 s”^. The two rate coefficients Ta^ and T^p 
are thus comparable to each other. 


V. THE INDIRECT REA IN LARGER 
MOLECULES 

For large polyatomic molecules, one expects that the 
electron capture into excited vibrational levels of the neg¬ 
ative ion M~ should be more efficient than for the CN” 
molecule. Nevertheless, the electronic non-adiabatic cou¬ 
plings along a single coordinate of M~ are not expected 
to be much different than in CN”. The latter statement 
should be true for molecules that do not exhibit singular 
effects near threshold, e.g. effects due to the presence 
of a resonance or virtual state, or if the radical ground 
state has an unusually large permanent dipole moment. 
In order to obtain another basis for comparison, we have 
also calculated the non-adiabatic couplings in the case 
of C 2 H- and found only slightly smaller values than in 







CN“. On the other hand, for molecules with large per¬ 
manent dipole (e.g. CgH, CgH or C 3 N), threshold effects 
could be important and the non-adiabatic couplings sig¬ 
nificantly larger. The study of the role of a permanent 
dipole moment in REA is out of the scope of the present 
study and has been discussed elsewhere 0 , 133 . 

Therefore, in the following development, we only con¬ 
sider the role of the number of degrees of freedom in the 
electron capture probability and discard the possible role 
of unusual threshold effects or vibrational Feshbach res¬ 
onances. Because the number of non-Born-Oppenheimer 
couplings and the energy-density Pm{Ec) of available vi¬ 
brational levels of M~ rises with the number of degrees of 
freedom, the overall capture rate is expected to increase 
correspondingly. One can readily estimate the density of 
final states in a molecule with several degrees of freedom. 
Let us thus consider a given linear molecular anion M~ 
and radical M with S degrees of freedom. Henceforth, 
the molecular ion M~ and its counterpart radical M will 
be represented as a set of vibrational harmonic oscillators 
with respective normal coordinates Qs and g', and asso¬ 
ciated harmonic vibrational frequencies Wg and w'. Both 
sets of normal coordinates represent the displacements 
around the equilibrium positions of the anion and neu¬ 
tral molecule, respectively. It will be shown later that 
for large enough carbon chain molecules, the harmonic 
oscillator model represents an accurate approximation. 
More importantly, we will not lose generality of the re¬ 
sults by using the harmonic oscillator approximation in 
our model. This is a good approximation near the min¬ 
ima of the M and M~ potential energy surfaces. 

In order to obtain a rough estimation of the density 
of vibrational levels, we first introduce a characteristic 
energy splitting AEc for the different vibrational levels. 
We then assume that the energy splitting in different 
modes qs is approximately equal to the average energy 
splitting, namely hiOg « AEc for all s (this approximation 
will be relaxed later). If we denote by Vg the number 
of excited quanta along the Qs coordinate, then the total 
number of excited quanta Vo corresponding to the energy 
of the system just above the radical threshold will be 
approximately given by: 


5 

Vo = Vs 

S=1 


Eea 
AEc ■ 


( 20 ) 


Appendix lA II) 


Nt 


V-H5- 1 
V 


( 21 ) 


Therefore, the density of vibrational levels is simply: 


Pm{Ec) 


1 /V -H 5 - 1\ 

AEc\ V J' 


( 22 ) 


The density of levels grows rapidly with the number of 
degrees of freedom and the total number of quanta. The 
largest contribution to the density comes from excitation 
of several modes at once, as shown in appendix lA 21 


A. Model to evaluate the IREA cross section 


Let us now introduce the model describing the IREA 
mechanism to large carbon chain radicals, neglecting the 
rotational motion. We need to evaluate the non-adiabatic 
couplings between the initial state of the M + e~ system 
and a vibronic state of M~ into which the electron is 
captured. It is convenient to introduce the dimensionless 
normal coordinates Qs = Qs \/Th0j7/^ of the negative ion 
M~, where Ps is the reduced mass. In a similar way, 
dimensionless coordinates qg' are also introduced for the 
neutral molecule M. For the molecules under consider¬ 
ation, the potential energy surfaces of M and M~ are 
nearly parallel, such that the normal coordinates qg and 
gs^ as well as the vibrational frequencies u>g and uj'g, are 
almost identical for all modes. This is the reason why we 
will use the same notations qg for the normal coordinates 
of M and M~. The only difference between the nor¬ 
mal coordinates of M and M~, which will be accounted 
for below in evaluating the Franck-Condon overlaps, is 
the displacement between the equilibrium positions of M 
and M~. Note that this approximation is not only con¬ 
venient, but has been shown numerically to be excellent 
for CN [l3], as well as for the hydrocarbon chains C 2 H, 
C 4 H and CeH 0. 

The non-Born-Oppenheimer operator acting on 

the nuclei coordinates and integrated over the electronic 
degrees of freedom, takes the form 


S 

A^’') = {■i>f\A\-i>ri\)r = '^Ag{qg) 

S = 1 


A 

dqg ’ 


(23) 


The affinity Eea of the interstellar anions ranges from 
about 3 — 4 eV, while the vibrational frequencies are 
usually of the order of 100 — 3000 cm“^. For the fre¬ 
quencies of the stretching modes, which are commonly 
the largest, this gives approximately 10 < Vo < 30. On 
the other hand, the frequencies of the bending modes 
are small {Hujg < 800 cm“^) and Vo could grow as large 
as Vo > 300 for very loose bending modes. For a given 
value of Vo, the number of vibrational levels of M~ in 
the interval of energy AEc is approximately given by (see 


where the electronic coupling Ag{qs) (equivalent to the 
operator introduced in Eq. (nm for the one-dimensional 
case) has the following value 

/ Q 

'V*f{r,^ — 'Vrix{r,^df. (24) 

We now assume that Ag{qg) are almost constant over 
the vibrational displacements, hence Ag{qg) Ri A^. This 
approximation has shown to be excellent in the case 
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of CN“ and C 2 H“, for which the electronic non-Born- 
Oppenheimer matrix element does not vary significantly 
over large inter-nuclei displacements (see Fig. El). We 
also assume that all the couplings have about the same 
value along the different modes, namely As{qs) ~ Ag for 
all s. With these approximations, the operator in Eq. 
(USD becomes simply 


A(’') 


5 


Ao^ 


A 

dqs ' 


(25) 


The calculation of the matrix elements of A^"^ requires 
the evaluation of overlaps between the initial and final 
vibrational wave functions, which cannot be expressed in 
analytical form in a general case. 

The matrix element of the non-Born-Oppenheimer op¬ 
erator (EH) to be evaluated is given by 

A/. = (x/A^Iw), (26) 

where \xi) = |00---0) denotes the vibrational ground 
state of M and \xf) = |uiU 2 "-U 5 ) denotes an excited 
vibrational resonant state of M~ . The vibrational wave 
functions are expressed as a product of harmonic func¬ 
tions of the qs displacements with Vs quanta in each nor¬ 
mal mode. If we now introduce the rescaled non-Born- 
Oppenheimer elements A/j = Afi/Ao, independent of the 
magnitude of the electronic non-adiabatic coupling, and 
insert Eq. (1251) into Eq. (1261) . we obtain 

A/, = ^(ui|0)..-(u,| —|0)...(U5|0). (27) 

s=l 


In the above sum, we call the mode s with the derivative 
d/dqs as the active mode, whereas all other modes will be 
referred as passive. In Eq. (EZl), the one-dimensional har¬ 
monic vibrational overlaps for passive modes take a sim¬ 
ple form. If we introduce the coefficients /3s = Aqs/V^, 
where A^s denotes the separation between the equilib¬ 
rium geometries of and M along the dimensionless 
normal coordinate qs , then the vibrational overlap is sim¬ 
ply given by [32| 


(^s| 0 ) = 



(28) 


We would like to stress the following points: 

• The value of /3s accounts for the frequency, re¬ 
duced mass, and length displacement of each mode. 
In terms of the original (not dimensionless) nor¬ 
mal c oordinate s Ags, these coefficients are /3s = 
Aqsx'p,sUJs/‘2h [Hj. 


• The sum over all vibrational quanta of the squared 
of the overlaps is unity, as required. 

The overlap, involving an active mode, is obtained by 
expressing the derivative operator using raising and low¬ 
ering operators 


In the case of a bending mode (/3s = 0) the expression 
(l29l) vanishes for all Vs except for Us = 1. It means that a 
degenerate active mode can capture an electron through 
non-Born-Oppenheimer coupling by an excitation of only 
a single quantum of vibrational excitation. 

The capture probability per unit time into the vibra¬ 
tional state \xf) of the anion is proportional to the square 
of the element of the non-Born-Oppenheimer operator in 
Eq. (Eni): 


lA/.P = 


S 

S = 1 


2vi 


•A" 


ui! • • • Us!■ 


■ - vs'. 


(/3f - 


2 e 


The coefficient /3o in the above equation represents the di¬ 
mensionless displacement between the minima of the neu¬ 
tral and anion potentials in the multi-dimensional space 
spanned by the S normal coordinates. It is simply given 

by /3o = \/^i H-+ /3? H-+4 • 

In our treatment, we did not calculate the positions 
and widths of the vibrational resonances, but instead 
computed the capture probability in each ” quasi-bound” 
vibrational state of M~ . Eor this reason, we now calcu¬ 
late an electron capture probability per unit time (P), 
averaged over a certain electronic energy interval Ia = 
[Etot — Al2,Etot -f A/2] around the total energy of the 
system. Then, the total capture probability per unit time 
is given by summing over all vibrational states whose en¬ 
ergies are situated within the interval Ia , namely 

4) = E lAl'- (30) 

/eiA 

where \Afi\^ is given by Eq. ([30|). In this approach, the 
probability in Eq. (I30p has to be evaluated numerically 
to obtain an exact value. However, (P) can be estimated 
analytically using the following arguments. Eirst, we fix 
the total number of quanta Vg introduced in Eq. (ISHI) . 
which is associated with an averaged total energy {E). 
Then, we evaluate the sum of transition probabilities to 
any vibrational states of AI~ with Vg excited quanta. 
Denoting Vf = vi + ■ ■ ■ + vs the total number of quanta 
of a vibrational state \xf) and introducing the sum 


• The bending modes have ds = 0 because no dis¬ 
placement exists between the minimum of the sur¬ 
face potentials of M~ and M along these modes. 
Therefore, the bending modes will only contribute 
negligibly in the electron capture. 


P(Vo)= ^ \Af,f, (31) 

Vj=Vo 

which is similar to the sum in Eq. (I30L but now evalu¬ 
ated fixing the number of quanta instead of the energy 
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interval, then the capture probability in Eq. (1301) can be 
approximated by the probability per unit time of exciting 
Vo quanta, namely 


In the above equation, lS.Ec is a characteristic energy 
spacing, taken as the difference between the average en¬ 
ergies corresponding to Vq and Vo + 1 excited quanta. 
The formula for 2i(Vo) is derived in Appendix iBl 


AVo) 


/3o^“”"(Vo+5-l)e-^2 

(Vo - 1)! 2 


(33) 


Choosing the vibrational ground state of M~ as the refer¬ 
ence energy, the average energy weighted over the capture 
probability to excite Vo > 1 quanta, is given by 


{E)^{Vo-l)h{uj), (34) 


where (w) is an averaged vibrational frequency (see Ap¬ 
pendix [B] for details). Therefore, the averaged cross sec¬ 
tion for electron capture in the IREA process of poly¬ 
atomic molecules is written in this model as 


Ac) 


A) 

Eel ^Ec 


AVo) 


TT Aq (Vq + S - 1) 2Vo- 2 -pl 

2Eei h{uj) (Vo - 1)! 


(35) 


where the current density in the incident plane wave of 
electrons, as well as the normalization of radial wave 
function, are the same as in Eq. m- 

The parameters entering Eq. (|55]) . namely the equi¬ 
librium positions, normal modes, average frequency 
and shift /Iq, should be calculated separately for each 
molecule. Below, we give typical values of these param¬ 
eters for the molecules under study. For bending modes 
in linear molecules, the shifts /3s between equilibrium po¬ 
sitions of M and M~ are equal to zero. Therefore, only 
stretching modes contribute into /3o. The average fre¬ 
quency of the stretching modes in large carbon chains 
is (w) ~ 2000 cm“^. Using a typical value of 4 eV for 
electronic affinity of carbon chain molecules C 2 nH and 
C 2 n-iN [n is an integer) and assuming a small energy 
of the incident electron, we obtain 13 < Vq < 16. Our 
calculated values of /3o are 0.19 (CN), 0.46 (C 2 H), 0.56 
(C 4 H) and 0.39 (CeH) [2^. The latter value, correspond¬ 
ing to the ground state of CeH, was calculated at 
the Hartree-Fock level, and could be somewhat underes¬ 
timated. In Fig. m the values of 2i(Vo) are plotted as 
a function of Vq for different values of /3o and choosing, 
for instance, the number of degrees of freedom in CgH, 
namely S — 22. As evident from the figure, 2i(Vo) is very 
sensitive to /3o but remains small for the typical /3o val¬ 
ues. As a result, the IREA cross section is small, unless 
/3o is much larger than unity, which is rather unlikely for 
other carbon chain anions that we have not studied. 



FIG. 8: Values of 2i(Vo) for 5 < Vo < 20 and S — 22. 


Taking the largest value of j3o = 0.56 obtained for 
C 4 H/C 4 H- with 5=10, (w) = 0.01 H, Vo = 15 and 
Aq = 0.13, which is similar to the one obtained for 
CN/CN“ (see Fig. [T]), we obtain 

/ N ... 5 X 10-17 2 

Ac) ~-Oq > (36) 

where the energy is in hartree (H). At electron en¬ 
ergy of 1 meV, the above formula gives the cross sec¬ 
tion of 10“ 1 ^ Oq. For somewhat more favorable situa¬ 
tion, when 5=22 (as for CgH) and a significantly larger 
/3o = 1 (an arbitrary value, not necessarily representing 
CgH/CsH-), the cross section at the same energy be¬ 
comes larger, about 10 “^ Oq, but still remains small to 
explain the astrophysical observations. 


B. Statistical distribution in the model 


We consider here the distribution of different anion vi¬ 
brational states after the electron capture if Vq is fixed. 
The electron capture probability into a vibrational state 
of the negative ion is proportional to the non-adiabatic 
coupling elements in Eq. pop . In this expression, 
j5l <C Vs in the factor (/3f — Vg) and, therefore, /3^ can be 
neglected. Making the change of variable —>■ Us -I-1 for 
the active mode in Eq. (1301) and summing over all the 
modes, we obtain 


|A/.p«(Vo + 5-l). 






ui! • • • ?;s! • 


■ • 


(37) 


with the new condition that ui -I- • • • -|- = Vo — 1- The 

most probable configuration {?;°, • • • vA to capture an 
electron corresponds to an extremum of the function 


S 

G{vi • • • ■us) = In \AfA - A(^ Us - Vo -f 1), (38) 

S = 1 
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where we use the method of Lagrange multipliers. Using 
the Stirling approximation, 

ln(us!) ~ Vs in i In(us) + In(v^), (39) 

the Lagrange equations, giving the extremum, become 

ln/32-ln?;° + ^-A = 0. (40) 


at least, for some molecules a rapid change of the vibra¬ 
tional state is possible. If it happens that the autodetach¬ 
ment width of this new state is small or the probability of 
the spontaneous emission for this state is enhanced, for 
example, due to a better Franck-Condon overlap, one can 
expect that the system will not be able to loose rapidly 
the attached electron and will eventually be stabilized by 
emitting a photon. If this happens, then the IREA cross 
section (HU) should be modified as 


Typically, /3f <C I and > 2, thus the third term in 
dMl) is much smaller than the first two terms and can, 
therefore, be discarded. Solving for we obtain 


V 


0 

s 


&(Vo-l) 

,"o 


(41) 


where the brackets denote the closest integer number 
because should be an integer number. The modes 
with the largest separation between the equilibrium po¬ 
sitions are, therefore, the most strongly excited in the 
electron capture, whereas bending modes barely con¬ 
tribute. However, if we consider the Sstr = N — 1 
stretching modes of a 7V-atoms carbon chain have simi¬ 
lar displacements 0s-, then the stretching modes will be 
equally populated with [(Vo — l)/5str] quanta. Since we 
are interested in large carbon chains Sstr > 4, the op¬ 
timum occupation numbers are such that Vg < 4. This 
result shows why the harmonic oscillator model is accu¬ 
rate for a carbon-chain molecule with several degrees of 
freedom. States with smallest possible number of excited 
quanta along each mode but involving many modes, con¬ 
tribute the most into the capture probability. Finally, 
only vibrational states with energies within the energy 
interval [(A) — AA; (E) + AE] with AE = HAuj^Vq — 1 
around energy {E) contribute significantly in the electron 
capture probability (see Appendix 1]). Because Aw ~ 
1000 cm“^, the energy spread is of the same order as the 
characteristic energy spacing AEc, such that the cross 
section of Eq. (1551) should give a correct estimate of the 
electron capture cross section. Nonetheless, because the 
affinity of a negative ion can usually not be assigned ex¬ 
actly to one Vo value, Eq. (l35l) only gives an order of 
magnitude of the cross section. 


C. Competition between radiative stabilization and 
autodetachment 

Eor large molecules, once the electron is captured, in 
addition to the processes of autodetachment and spon¬ 
taneous emission from the formed resonance, the system 
can also change its vibrational state. Although the en¬ 
ergy of the system cannot change without emitting a 
photon, such a change in the vibrational excitation is 
possible because all vibrational levels in this energy re¬ 
gion are in fact resonances and have finite widths. In 
large polyatomic molecules the coupling between vibra¬ 
tional modes could be strong, so one would expect that. 


criREA = --CTc , (42) 

tot 

where Fjot now includes also the width r„c for the men¬ 
tioned change of the vibrational state F tot = Lsp -I- Fad -I- 
r„c- The time scale for change of vibrational state in 
a favorable case when the initial and final vibrational 
states strongly interact with each other, must be of the 
order of a vibrational period, in a ps range, which makes 
r„c dominant over the two other widths. In this limiting 
case, the IREA cross section (Jirea would be determined 
by its upper bound CTc of Eq. dni). Based on the above 
discussion, we can conclude that in the process of for¬ 
mation of negative ions {m = 1,3,5) and C„H“ 

(n = 2,4) by radiative attachment, a significant contri¬ 
bution into the total REA cross section of the IREA is 
unlikely. 

Non-adiabatic couplings in polyatomic molecules could 
be significant near a conical intersection of potential en¬ 
ergy surfaces and, as a result, lead to a larger probability 
of electron capture during the first step of IREA. Eor the 
considered case of radiative attachment, this can happen 
for certain incident electron energies if (a) the neutral 
molecule has degenerate vibrational modes and (b) the 
formed electronic state of the anion is also degenerate. 
The linear carbon chain molecules CnH {n = 2,4,6,8) 
and CmH (m = 3,5) do have degenerate vibrational 
modes, and some of the corresponding anions may have 
excited electronic states of the degenerate ^H symmetry. 
If such excited electronic states have appropriate ener¬ 
gies, an electron with tt symmetry, incident at the neu¬ 
tral target, could be captured due non-adiabatic Renner- 
Teller coupling, into the degenerate ^H state of the an¬ 
ion, exciting at the same time, the degenerate vibra¬ 
tional mode of the molecule. Such a process is similar to 
electron capture in dissociative recombination of linear 
molecules, HCO+ and N 2 H+ [33l - [^ . but there is an im¬ 
portant difference: In collisions between an electron and 
a positive molecular ion, the density of electronic reso¬ 
nances is larger by several orders of magnitude. In dis¬ 
sociative recombination, electronic Rydberg resonances 
with a vibrationally excited core could be found virtu¬ 
ally at any energy above the lowest ionization limit. But 
only a few molecular anions, such as C^, are known to 
have electronic resonances. Therefore, for a significant 
increase of the IREA cross section, the anion should have 
an electronic ^H state close to the energy of the initial 
continuum state of the e~ — M system. This could, prob¬ 
ably, be the case for one of the carbon chain anions, but 
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it seems unlikely that all observed anions have an elec- this manuscript is based on work conducted while A. Orel 
tronic resonance near the energy of the initial state of the was serving at NSF. 
e~ —M system with small energy of the incident electron. 


VI. CONCLUSION 

We have extended the theoretical approach to study 
the process of radiative electron attachment, developed 
in our previous study [s^, to larger molecules. Using the 
approach we have calculated REA cross sections for the 
three negative molecular ions CN“, C 2 H“, and C 4 H“. 
The following concluding remarks should be stressed as 
a result of the present study. 

• For completeness of the approach, two pathways for 
the process have been considered: (1) In the direct 
pathway, the electron, incident on the molecule, 
emits a photon and forms a negative ion in one 
of the lowest vibrational levels. (2) The indirect 
pathway is a two-step process, for which the in¬ 
cident electron is initially captured through non- 
Born-Oppenheimer coupling into a vibrationally 
excited state of the anion, forming a resonance. As 
a second step in IREA, the resonant vibronic state 
of the anion emits a photon, which stabilizes the 
anion with respect to autodetachment. The con¬ 
tribution of the indirect pathway was found to be 
negligible compared to the direct mechanism if no 
unusual threshold effects, virtual states or vibra¬ 
tional Feshbach dipole resonances are present. 

• The obtained REA rate coefficients evaluated at 
temperature 30 K are 7 x 10“^®cm^/s for CN“, 
7x 10-i^cm3/s for CzH”, 2 x IQ-i^cm^/s for C 4 H-. 
The coefficients depend weakly on temperature be¬ 
tween 10 K and 100 K and increase relatively fast 
with temperature above 200 K. The validity of 
the obtained results is verified by comparing the 
present theoretical results with experimental data 
from recent photodetachment experiments. 

• Previously, it was believed that carbon-chain an¬ 
ions, C„H“ {n = 2,4,6, 8 ) or (m = 3,5), 

observed in the interstellar medium, are formed by 
the REA process. The REA rate coefficients ob¬ 
tained in this study are too small to explain the 
observed abundance of the anions in the ISM. For 
example, for C 4 H“, the magnitude of the rate co¬ 
efficient needed to explain the observed abundance 
should be of the order of 10“^° cm^/s Q. Thus, 
the present results suggest that in the ISM, either 
dipole resonant states or non-local threshold effects 
increase drastically the REA cross section (and the 
rate of anion formation) or these anions are formed 
through a process different than REA. 
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Appendix A: Formula for density of states 


1. Total number of states 


We assume that we have S vibrational modes of the 
negative ion M~, and the energy of quanta in different 
modes is approximately the same AE. To reach energy 
Ei + Eel of the initial state of the M + e~ system, an 
excitation of V quanta is needed 

S 

Vs =V with Vs >0 . (Al) 

S=1 


The density of vibrational states of the M~ anion near 
energy Ei + Eei is evaluated as Nt/AE, where Nt is 
the number of combinations {ui, U 2 , • • • , I’s} how the sum 
above can be formed. 

To obtain the number of combinations Nt, it is conve¬ 
nient to represent V quanta as objects arranged in a row. 
The row of quanta is separated in subsets vi,V 2 , - ■ ■ ,vs 
by iS — 1 walls. Now we would like to count the num¬ 
ber of possibilities how the 5 — 1 walls can be placed. 
We can place the first wall in V -I- 1 places between the 
objects, keeping in mind that we can also place it be¬ 
fore the first or after the last quantum. The second wall 
can be placed at V -I- 2 different locations, because now 
we have V quanta plus the first wall in the row. (We 
will account later for the fact that the walls are identi¬ 
cal) . Continuing in this way, we obtain the total number 
S-i 

n (V + f) = (V -I- 5 — 1)!/V! of possibilities of how the 

i—1 

5—1 distinguishable walls can be placed. Since the walls 
are all the same, we have to divide the product by the 
number (5 — 1)! of permutations of the (5—1) walls for 
a given partition of the row of V. Thus, the number of 
combinations corresponding to the sum of Eq. (ED is 


N'y' 






(A2) 


2. Number of density of states when only n modes 
are excited 

We derive now the formula for the number of combina¬ 
tions corresponding to n excited modes, given that there 
are 5 modes. For example, if only the first n modes are 
excited, we have 


n 

Vs =V with Vs > 0. (A3) 

S=1 

In this situation, the number of combinations how the 
sum above can be formed is given by ■ To prove 
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this it is again convenient to represent the V quanta as 
objects arranged in a row. The row of quanta is separated 
in groups of ui, U 2 , • • •, Vn objects, by n — 1 walls. Now 
we would like to count the number of possibilities such 
that n — 1 walls can be placed. We can place the first wall 
in V — 1 places between the objects, keeping in mind that 
at the left and the right of the wall there were at least 
one quantum because Vs > 0. The second wall cannot 
be placed at same the position, because it would mean 
that one of the n modes will be inactive, i.e. Vs = 0 . 
Therefore, for the second wall there are V — 2 ways to 
place it. The total number of possibilities to place n—1 

n— 1 

walls is, therefore, n(v — i). Since the walls are all 

2—1 

the same, we have to divide the product by the number 
(n — 1 )! of permutations of the identical (n— 1 ) walls 
for a given partition of the row of V objects. Thus, the 
number of combinations corresponding to the sum of Eq. 
(Tnil is indeed 

To obtain the total number of combinations Nn when 
any n modes are excited, we have to multiply the above 
number with the number of ways how n modes can be 
chosen from the set of S modes. There are (^) such ways. 
Therefore, 


Nn = 



(A4) 


Appendix B: Closed-form expression for I(Vo) and 
statistical formula 


We will reduce T(Vo) in Eq. (|3T]l to a closed-form ex¬ 
pression. In order to simplify the equations, we introduce 
the shorthand notation Iq = 2I(Vo) exp(/3Q). Erom Eq. 
O, we decompose Iq in three terms: 


io 


E E 


s=l V/=Vo 






/HD- 


(Bl) 


The under braces indicate the factors involved in each 
term, such that Iq = Ti -f 2 I 2 + Ta ■ For convenience, we 
introduce the function 


Z(/3i--./35)= ^ 

V/=Vo 


p2vi...p2vs...p2vs 

Vi \ ■ ■ ■ Vsl ■ ■ ■ Vsl 


(B2) 


Recalling the multinomial theorem, this function takes 
the compact form 


Z{Pi---Ps) 


(/?? + ••• + 

Vo! 


Vo! ■ 


(B3) 


The first term Xi is expressed as 


Finally, we will prove that the sum of Nn gives the 
total number Nt of vibrational states of Eq. (IA2I) . i.e. 
we will show that 

S 

NT = Y,Nn. (AS) 

n—1 

It can be verified using Vandermonde’s identity, which 
states that 


5 

ii=^^2^(/3i...^5), (B4) 

S=1 

such that we readily obtain 

o2Vo+2 

The second term I 2 is written as 



In our case, we have 



where S' = 5 — 1. Now using Vandermonde’s identity 
with TO = V— 1, /=5'-|-l and r = 5', the sum is 
written 




n—1 


fv + s-i\ 

I 5-1 


(AS) 


T 2 


E E 

s=l V/=Vo 


Ui! • • • Us! • • • U5! 


(B 6 ) 


which takes the convenient form 

j _ _ ■ ■ ■ Ps) 

^ ^ 5 In /3s 


(B7) 


Using the value of Z(/3i • • • Ps) in Eq. (|B3I) . we obtain 
the following expression: 


12 


2 Vo^g^° 

Vo! 


(B 8 ) 


The last term I 3 is clearly the dominant one, with the 
slightly more complicated form: 


is 


5 

Ui!---Us!---U5! 


which is indeed Nt of Eq. (|A2p . 


(B9) 
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Note that the apparent indeterminacy is easily removed 
for a degenerate mode by taking the limit /3s —?► 0 in Eq. 
(IB9() . Applying a similar method, we find 


- _ 1 \ ^ 1 d'^Z{Pi ■ ■ ■ Ps) 
^ 4^/32 i9(ln/3s)2 


(BIO) 


and thus we conclude that 

- Vo (Vo+ 5 

4-3 =-ry-i- 


(Bll) 


Finally, the closed-form formula for 2i(Vo) is given by: 


I(Vo) = 


[(Vo - /3g)^ + Vo(5 - 1)] 
Vo! 



In the above expression, Pq Vo, which allows us to 
simplify the expression as 


I(Vo) 


(Vo + 5 - 1)/32^‘>-2e-^o 

(Vo-1)1 2 


(B13) 


We also derive an expression for the average total en¬ 
ergy of the system, weighted over the electron capture 
transition probabilities, given that Vo quanta are excited. 
First, neglecting the terms Pg <C Vg in Eq. (|5(I1) . the 
transition probabilities to capture in an anion vibrational 
state are proportional to 


r{Pi---ps) 


p2vi...p2v.-2...p2vs 
Ui! • • • (Ws - 1)1 • • • 


(B14) 


Therefore, the partition function used to normalize the 
probabilities takes the form 

5 

S = ^ (B15) 

s=l V/=Vo 

and the average energy is thus expressed as 

(^) = 4 E EfV{P,---Ps), (B16) 

“ Vf=Vo 


where Ef = h{viUJi -I- • • • -I- vsojs) is the energy of the 
vibrational state \xf) (the anion ground state is chosen 
as the reference energy). After some straightforward ma¬ 
nipulations, we obtain for Vo > 1 the expression 

(E) ^ iVo-l)h{uj), (B17) 

where (uj) = V) Ps^s is an average vibrational fre- 
quency weighted over the displacements /3s. In a simi¬ 
lar fashion, one can calculate the energy spread AE = 
yj (F^) — (i3)2 around the average energy (A), which 
takes the value 


AE Ri hAuj^/Vo - 1, (B18) 


where Aw = Y^(aVy^^'(wp^ is the vibrational frequency 
spread, with (w^) = ^Ws- 
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